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^^ We derive an exact formula for the complex frequency in spatio-temporal stabil- 

^^ ity analysis that is valid for arbitrary complex wave numbers. The usefulness of the 

&J0 formula lies in the fact that it depends only on purely temporal quantities, which 

^ are easily calculated. We apply the formula to two model dispersion relations: the 

■^^i linearized complex Ginzburg-Landau equation, and a model of wake instability. In 

f~^ the first case, a quadratic truncation of the exact formula applies; in the second, the 

T-H same quadratic truncation yields an estimate of the parameter values at which the 

transition to absolute instability occurs; the error in the estimate decreases upon 

;-H increasing the order of the truncation. We outline ways in which the formula can be 

^^ used to characterize stability results obtained from purely numerical calculations, 

'"O and point to a further application in global stability analyses. 

I 
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^ 1. Introduction 

C/3 We are concerned here with linear stability analyses of general fluid flows, and 

^^ begin by recalling the main steps of these, before introducing the motivation of 

the present work, thereby also introducing the necessary notation. We shall first 
restrict ourselves to the linear stability analysis of a parallel flow Uq{z) to two- 
dimensional disturbances; the extensions of the work to non-parallel flows and to 
three-dimensional disturbances are set out in Sec. [6} Such analysis is accomplished 






> 

^N^ tion at fixed parameter values (such as the Reynolds number Re) yields an eigen- 

t^ — value problem that connects the wave number, the frequency, and the growth rate 



by solving the Orr-Sommerfeld (OS) equation (fOrr 1907). The solution of the equa- 



^~~l of the disturbance ( Drazin and Reid[ 1981 Orzag 1971). For two-dimensional dis- 

(T^ turbances, the perturbation streamfunction has the form 'ip{x, z, t) = (/)(z)e'("^~"*\ 

^^ where a = a-c + iai and uj = uj-c + iwi are complex numbers, a-c is the (real) wave 

Cn number, lo^ is the (real) frequency, and ai and Wj are the spatial and temporal 

. . growth rates respectively; in this context, the OS equation reads 

> 2 

•'-j \a{Uo{z)~c){dl-a^)(t)-iaUl!,{z)(j)^Re-^{dl-a'^) <p, c = w/a. (1.1) 

5h That the wavenumber a is complex indicates that the instability is spatio-temporal 



in nature (Huerre and Monkewitz 1990); in a temporal analysis (Drazin and Reid 
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1981), tti = and the eigenvalue problem is solved for aj; and uj-c. This results in 



the following temporal dispersion relations for the temporal frequency and growth 
rate respectively: w'^^p = a;r(ar,ai = 0), and uj-^™"^ = a;i(Q;i.,ai = 0). The basic 
flow is called unstable if the growth rate is positive for some ai-. It is of interest 
further to classify unstable flows as convectively or absolutely unstable: the flow is 
convectively unstable if initially localized pulses are amplified in at least one moving 
frame of reference but are damped in the laboratory frame; on the other hand, the 
flow is absolutely unstable if such pulses lead to growing disturbances in the entire 
domain in the laboratory frame. 



Briggs (1964) developed a method to classify unstable plasmas according to this 



dichotomy. This criterion was extended to fluid dynamics by |Huerre and Monkewitz] 
(1990) (see also Lingwood (1997); Schmid and Henningson (2001)). The approach 



is based on an eigenvalue analysis of the OS equation: (i) a necessary condition for 
absolute instability is if the imaginary part of the frequency is positive at a saddle 
point in the complex a-plane: Wjo '.— a;i(ao) > 0, where ag solves du/da — 0; 
(ii) to obtain a sufficient condition for absolute instability, the saddle point ao 
in the complex a-plane must be the result of the coalescence of spatial branches 
that originate from opposite half-planes at a larger and positive value of Wi (this 
coalescence or 'pinching' of the spatial branches is accompanied by the formation 
of a cusp at Wjo in the complex lo plane) . Typically, these steps can be implemented 



efficiently using numerical spectral methods (Boyd 2001). 



An analytical connection between spatio-temporal and temporal growth rates 
is desirable for a number of reasons . First , in flows where several temporal modes 
are unstable (e.g. O Naraigh et al. (2011)), the interpretation of the results of a 



spatio-temporal stability analysis becomes difficult, as does the mere determina- 



tion of the convective/ absolute (C/ A) transitions ( Suslov , 2006 ) . In this case, w- 
plots of the distinct spatio-temporal modes consist of complicated interpenetrating 
surfaces, whose domain is the entire complex a-plane. On the other hand, many 
theoretical tools exist for characterizing temporal instabilities, including energy- 



budget analyses (e.g. Boomkamp and Miesen (19961) and asymptotic analytical so 



lutions (Drazin and Reid 1981), so the governing physical mechanisms for growth 



of temporal modes is often well documented and competing temporal modes can 
easily be distinguished. A direct link with temporal mode could therefore enable 
one to distinguish the competing spatio-temporal modes. Furthermore, in para- 
metric studies involving large parameter spaces (e.g., thermal boundary layers or 
two-phase flows), wherein the C/A transition curves are plotted as a function of the 
flow parameters, the dependence of the C/A transition curves may be difficult to 
analyse. For example, for a system with a two-dimensional parameter space (/^i, /i2), 
the curve demarcating the transition between convective and absolute instability is 
given by the generic formula /i2 = fifJ-i)- Here, knowledge of the physical proper- 
ties of the system (in particular the familiar temporal stability properties) may be 
beneficial (if not crucial) to deduce the function /(•). 

The formula proposed herein satisfies these criteria. Gaster ( |1962 ) made some 
progress in this area for purely spatial instabilities. Such an instability corresponds 
to setting Wj = in the OS eigenvalue analysis, and computing the resulting spatial 
growth rate ai as a function of (wr,ar). Practically, this corresponds to a localized 
disturbance in the flow that oscillates at a characteristic frequency cj,., and grows 
downstream of the disturbance at a (spatial) rate —ai. Spatial growth therefore cor- 
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responds to negative values of ai. Using the fact that the general complex frequency 
Lo{a) is an analytic function, Gaster integrated the resulting Cauchy-Riemann con- 
ditions along a straight-line path of small length, (ar,0) — >■ (ar,ai), with |ai| <S^ 1. 
The resulting integrands are thus regarded as constant and, setting Wi — 0, Gaster 
obtained the following formula for the spatial growth rate: 

temp /duj'^'^P 

This equation provides a precise connection between the spatial growth rate and 
more familiar temporal quantities. Nevertheless, it is limited, in the sense that 
it applies only to spatial growth, and is valid only for small values of |ai|. The 
specific aim of this paper is to start with the Cauchy-Riemann conditions for uj{a) 
and obtain an exact formula that connects the spatio-temporal growth rate to 



purely temporal quantities. Gaster's formula (1.2) emerges from this approach as 
one particular limiting case. 

This work is organized as follows. We derive an exact formula connecting the 
spatio-temporal growth rate with purely temporal quantities in Sec. [2J We dis- 
cuss the two lowest-order truncations of this formula in Sec. |3j Gaster's formula, 
and a situation in which the temporal quantities depend only quadratically on the 
wavenumber ctr. This enables the formulation of condition (i) into a succinct equa- 
tion encoding the competition between in-situ growth and convective effects. In 
Sec. [4] we describe the singularities that typically occur in the dispersion relation, 
and discuss how these can hamper the convergence of our formula. In Sec. [5] we 
apply our results to simple test cases and describe further applications of the exact 
formula and its quadratic approximation. 

2. An exact formula for the complex frequency 

The development of the formula starts with the assumption that an eigenvalue anal- 
ysis of the OS equation yields a complex frequency uj that depends on the streamwise 
wave number a as an analytic (holomorphic) function, uj = a;r(ar, Q^i) + iwi(ar, ai), 
where we have fixed the other system parameters (e.g. Reynolds number). The do- 
main D of analyticity is an open simply-connected subset of Z? C C, containing a 
part of the real line. Details concerning the domain are discussed further below. As 
a consequence of the analyticity of uj{a), we have the following Cauchy-Riemann 
conditions: 

;^— = ^— , (2.1a) ^~ = ^^~- (2.1b) 



We now use Eqs. (2.1a| and (2.1b) to derive an expression for Wi(ar, cki) as a function 



(2.2) 



of the purely temporal stability properties. 

We work at a fixed value of a-n in Eq. (2.1a) and identify the group velocity 
Cg(ar) = (duj^/dar)(a^fl)- Furthermore, we Taylor-expand duj^/da^ into the region 
D: 

a^ = 2. ^^s«K)ai , c,„(a.) = — — 

n— ^ ai— u 

This amounts to a complex Taylor expansion centred at (ar,0), and is therefore 
valid on a disc of radius R contained entirely in D. The magnitude of the radius of 



Article submitted to Royal Society 



L. O Ndraigh and P. D. M. Spelt 



convergence R is discussed below. Each term in this Taylor expansion is available 
from a purely temporal analysis. We describe this process for the first few terms 
(for brevity, we use di := d/dai and d^ := d/da,-): 

1. n = : Cgo = 9iWr|Qi=o — Cg{a^), the group velocity. 

2. n = I : We have the following sequence of steps: 

72 temp 
Ccri = didrOJr = d,:diljJr = ' -^i^rWi "= -^ , (2.3) 

where Wj"™^ is the purely temporal growth rate. 

3. n — 2 : We have 



da'^ 

(2.4) 



4. n = 3 : We have 



Cg3 = didid-Auj^ = didididiUJr = ' —dididrdrUJi = —didrdrdiUJi 

j4 temp 

= ' —didrdi-drUJi- = —d^dydrdiUJr =' di-di-drdrCjJi — — r-j — -(2.5) 

daf 

5. n — 4 : We have 

Cg4 = dididididiLUi- ~ dididid^diuj^ = ' —didididi-d^uji ~ —d\d\d^d^diduji 

«J- a,.a,a,.9,9,..,"^°^. (2.6) 

We deduce the higher-order terms by inspection. Assembling these results, we obtain 
an expression for duj-c/da^, extended into the complex plane: 






K,oO ;^o (2n)! da?""' +Z.(2n + 1)! da?»+2 «i • ^ ■') 



Next, we use the Cauchy-Riemann condition (2.1a) and connect the Taylor ex- 



pansion in Eq. (2.7| to an expression for Wi(ar,ai), valid in the entire complex 



half-plane a-c > Q. We have 



n J2n„ °° ( T\n+1 j2ri+2, temp 



LjJ: 



aai aa,. ^ (2n)! da?" ' ^^(2^ + 1)! da?"+2 ' ' ^ ' 

n— n— 

We integrate this relation along a line perpendicular to the real axis, from (a,., 0) — > 
(ar,ai) (ai is arbitrary, i.e. |ai| is not necessarily small). The result is 

°° ( ^\n j2n„ °° ( i Vn+l j2n+2, ,tomp 

I \ tomp/ N , Y^ \ ^1 " Cg 2n+l , V^ I J-J " "-'i 2n-|-2 

.K«r,aO=-, K) + L(2,+ l),d«2„"i +2.(2^ + 2)! da?"+2 '^^ ' 

n— n— 

(2.9) 
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This is an exact formula for the imaginary part of the complex frequency, which 
depends only on purely temporal quantities. Mathematically, this is a trivial state- 
ment: since the complex frequency is analytic, analytic continuation implies that its 
behaviour on the real line completely determines its behaviour on D. Nevertheless, 
this formula will help us to create a practical means of predicting the threshold for 
absolute instability. 



The series (2.9) amounts to a complex Taylor series centred at the point (ar, 0) 
and therefore converges inside a disc of radius R. The radius R is the minimum 
distance from the point (ckr,0) to the nearest singularity of io{a). Care must be 
taken that the point of interest (e.g. the saddle point of (^{a)) lies inside this disc. 
A detailed discussion of the radius of convergence of Eq. (2.9 1 is given in Sec. |4] 



Low-order truncations and an approximate criterion for 

the C/A transition 



In this section we investigate the implications of truncating Eq. ( 2.9 1 at linear order 



and at quadratic order. At linear order, the spatial analysis of Gaster is recovered; 
at quadratic order, the complex Ginzburg-Landau model is obtained, albeit with 
further understanding, in the form of a 'balance condition' for the onset of absolute 
instability. 



(a) Relation to the analysis of Gaster 



The formula of Gaster ( |1962 ) can be recovered directly from the approach in 
Sec. [2] Consider Eq. (2.9) with |ai| <^ 1, such that only the lowest-order term gives 
a contribution to the equation. The result is 



ujiiar, ai) ~ Wi™P(ar) + Cg(ai.)ai, 



"i 



0. 



(3.1) 



On a spatial branch, we have Wj = 0, hence Eq. (3.1 ) becomes a; ~ —LL)-^'^"^^{a^)/cg{a^) 



as |ai| -^ 0, which is precisely the formula of Gaster (1962) (Ec^. ( |1.2[ )) for spatial 
modes. 



(6) The quadratic approximation 
We also derive conditions for the C/A transition based on a second-order trun- 



truncation, the Taylor series (2.9) reduces to 



cation of the series dispersion relation (2.9) - the quadratic approximation. In this 



j2 temp 
_ temp/ N I „ /„, \„. 1 " W; ^2 



Wi(ar,ai)=H P(ai.) + Cg(ar)ai 



da? 



(3.2) 



The necessary saddle-point condition for an absolute instability (as reviewed in 
Sec. [1]) is duji/da^ — dijj\/da\ = 0; using Eq. (3.2), these conditions amount to 



dc^,'°"^P dc. , d3c^,"='"P , „ . . d2 temp 



-tti 



dofr ' da^ ' ^ daj? ' 
Taking the second condition, we get 



0, 



Cg(a,.) 



dai 



-ai = 0. 



Cg(ar) 



dai 



(3.3) 



(3.4) 
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Substitution into the first condition yields 



r{a. 



dcg, 
dar 



dc^f'""P ^2 temp 



dar 



da? 






K) 



d3 



W: 



temp 



da? 



'dV^ 
da'? 



(3.5) 



(Note that third-order derivatives appear in this calculation - but only through the 
process of finding the saddle-point location, and not in the series expansion of the 
dispersion relation.) 



The existence of a real root of Eq. (3.5) is a necessary condition for a saddle 



point to occur. Once the root of Eq. (3.5) is extracted, we derive a condition for 



the Wi to vanish at the saddle point, as this is the sign of the transition to absolute 



instability. Referring to Eq. (3.2), uj[ vanishes at the saddle when 



d2w,"='"P 



da'? 



w*""P(a: 



2^g 



(«:)• 



(3.6) 



where a* is the root of Eq. (3.5). Knowledge of the saddle-point location a*, to- 



gether with the dispersion relation (3.2) yield further information that can be used 



to verify if the saddle point arises as a result of the coalescence of spatial branches 
ramifying into different halves of the complex a-plane (i.e. the necessary and suffi- 
cient conditions for the onset of absolute instability) . 

The dispersion relation for the linearized complex Ginzburg-Landau equation 
is quadratic in the wavenumber, and no error is incurred in this case in applying 



the quadratic truncation of the series dispersion relation (2.9) 



cj = Ua + ^Coa^ 



it^l 



a e 



(3.7) 



here f/, cq, and /xq are positive real constants. Using the direct approach (Huerre 



2000), the system (3.7) is absolutely unstable when fiQ exceeds a threshold value: 



where the saddle point occurs at wavenumber 

a = -[/ (co + i) / (1 + eg) . 



(3.8) 



(3.9) 



We compare these results with the framework developed here in this section. The 
group velocity for the dispersion relation (3.7) is Cg = U -\-coar, while the temporal 
growth rate is ui^'^™^ — (/Iq — a^)/2. Substitution of these two formulae into the 



saddle-point condition (3.5) yields a value for a^- at the saddle in exact agreement 



with Eq. (3.9). Substitution of the same ar-value into condition (3.6) yields the 
relation /ig = U{1 + Cq)~'^ for the onset of absolute instability, in exact agreement 



with the known condition (3.8) 



An advantage of the framework developed in this section (compared with the 



direct approach) concerns the insight given by Eq. (3.6) into the balance of compet 



ing effects that brings about absolute instability. Intuitively, we may think of the 



condition (3.6) as a competition between in-situ growth on the left and convection 



of the disturbance on the right. Absolute instability can set in only if the in-situ 
growth equals or exceeds the tendency of the parallel flow to convect disturbances 
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downstream. Finally, we note that although a 'quadratic approximation' was pur- 
sued in earlier work (Deissler 1987 Monkewitzj [1988 ), this was in relation to a 
numerical, iterative procedure for determining the location of the saddle point, and 
did not contain an explicit condition linking the C/A transition to the temporal 
analysis. 



4. Convergence of series for cui in the presence of a bestiary 

of singularities 

In this section we review the different kinds of singularity that occur in the disper- 
sion relation w = u}{a), and discuss the implications of this non-analytic behaviour 



for the convergence of our formula (2.9). Here, the 'impulse response function' refers 



to the solution of the spatiotemporal problem 
[Uo{z)d, + dt] (dj+dl) i^{x,z,t) - Ui;{z)d^^{x,z,t) = i?e-i (9^ 

subject to the in initial condition 

ip{x, z,t — 0) = S{z)S{x). 



dl) '4){x,z,t), 



(4.1) 



Equation (4.1) is typically solved using a Laplace-Fourier decomposition; the dis- 



persion relation of the resulting eigenvalue problem (the OS equation) determines 
the shape of the impulse response. 

1. Logarithmic branch cuts associated with a critical layer: In inviscid flows, the 
critical layer z^ for which t/o(^c) — c induces a logarithmic singularity in the 



dispersion relation; see Huerre and Monkewitz ( 1985 ) . This in turn splits the 



impulse response function into two parts: a continuous part, associated with 
the logarithmic singularity, and a discrete part, associated with the zeros of 
the dispersion relation. The continuous part is not important, as it decays al- 



gebraically in time (Huerre and Monkewitz 1985 1; thus, the focus remains on 



the discrete part. This singularity anyway disappears for the finite- Reynolds- 
number case, wherein the coefficient of the highest-order derivative in the OS 
equation never vanishes. 

Branch cuts in unconfined flows: In so-called unconfined flows, the disper- 
sion relation w = uj(a) may possesses branch cuts along the imaginary axis 
Qfr = ( Huerre and Monkewitz[ 1985[ Lingwood 1997). In such systems (e.g. 
mixing layers, flow past obstacles), both the streamwise and normal dimen- 
sions of the fluid container extend to infinity (labelled by coordinates x and z 
respectively), and the disturbance streamfunction decays to zero as \z\ — > oo. 
Then, the asymptotic (|2;| — )■ oo) OS equation reads 



ia (f/oo - c) (af -a^)(t)^ Re-^ {dl - a^) cj), 



:a;/a. 



for a symmetric base state Uo{z) with 



Uo, 



lini Uo{z)^ and 

:— >±00 



hm ui;{z) 

:— >-±oo 



This asymptotic solution possesses an inviscid mode 

g-sign(7,)7^^ where7= ^/, 



g--sign(a,)az and 



viscous mode 



'a^ 



iaRe{Uoo — c). The inviscid 
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mode induces a branch cut along the imaginary axis (Huerre and Monke 



witz 



1985) 



while the viscous mode induces hyperbolic branch cuts in the 
a-plane ([Aships and Reshotko 1990 Schmid and Henningson 2001[). All of 



these branch cuts (or 'continuous spectra' ( Schmid and Henningson } 2001)) 
contribute to the contour integral associated with the impulse-response func- 
tion. The contour integral of the impulse-response function therefore possesses 
two contributions: one from the zeros of the dispersion relation uj{q) = 
('poles'), and the other from the continuous spectrum. However, to diagnose 
absolute instability, it suffices to consider the discrete part, as the continu- 
ous spectrum can produce temporal growth only when the discrete part is 



absolutely unstable (Lingwood 1997) 



3. Discrete poles along the imaginary axis in the complex a-plane: Healey (2007 



2009 ) has shown that confining an inviscid version of the flow described in (2) 



between two plates at z — ±H causes the character of the singularity along 
the axis ar = to switch from a continuous branch cut to a set of discrete 
poles. For, the inviscid OS (Rayleigh) equation 



; {Un{z) - c) (92 - q2) - iaUo{z)<j) = 0, 



can be re-written as 



{dl 



a 



ia^o(^)'/' 
iaUo{z) — iw' 



such that the singularity a; = 00 corresponds to {d^ — a^)0 = 0. This lim- 
iting equation satisfies the appropriate boundary conditions for a confined 
fiow (no-penetration at the walls where confinement is enforced) whenever </> 
can be made into an oscillatory function, that is, when a = -imiri/ H^ with 
n = 1, 2, • • • . These poles induce a series of saddle points near the imaginary 
axis. These saddle points can dominate over the saddle point associated with 
the analogous unconfined flow for sufficiently small H . Such saddles can even 
pinch (according to the Briggs criterion) and thus produce absolute insta- 



bility (Healey 2007 2009), leading to the conclusion that a confined flow is 



'more absolutely unstable' than its unconfined analogue. 

Since the introduction of viscosity to the problem does not change the large-w 
form of the eigenvalue problem, it is expected that such poles will persist in 
the viscous case (with some modification of the large-w streamfunction arising 



from the different boundary conditions in the viscous case) (Healey 2012). 
This contention is supported by numerical evidence (Sec. pi). 

Further branch cuts in the complex a-plane: The generic dispersion relation 
2?(a,aj) = of the OS equation naturally evokes consideration of the mul- 
tivariable function I?(a,a;) itself. This is an analytic function in each of its 



variables (Caster 1968). Moreover, dV/du has at least one zero (Gaster 



1968). In the neighbourhood of this point (labelled (ao,wo))i 2?(a, w) can be 



expanded as 

'D{a,uj) = 'D{ao,ujo) +T)a{ao,ujo){a - ao) -\- ^V^ 



+ Higher-order terms. 
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where the subscripts a and w indicate partial differentiation. If the points 
(a,a;) and (aojWo) satisfy the dispersion relation 2? = 0, then this equation 
can be recast as 



uj — ujo + 






1/2 



{a - ao) 



1/2 



(4.2) 



The dispersion relation uj{a) is therefore non-differentiable at ao and possesses 
a branch cut along a line emanating from the point ao- Equivalently, one may 
regard the dispersion relation as being multi- valued, taking distinct values on 



the two distinct sheets of the Riemann surface described by Eq. (4.2) 



If a saddle point of the dispersion relation is sufficiently close the branch 



cut (4.2 1, the branch cut can prevent the saddle point from pinching, such that 
the Briggs criterion is not satisfied (see SecfTl). Equivalently, in such cases, it 
is not possible to deform the contour in the impulse-response integral into a 
steepest-descent path without enclosing the singular point ao. Such saddles 
do not therefore generate absolute instability (Lingwood 1997). 



It is clear from this list that knowledge of the global topography of the dispersion 



relation is required to determine absolute instability (on this point, see also Ling- 
wood (1997)). Equally, such knowledge is necessary to determine when Eq. (2.9) 
can be used in confidence. Since Eq. ( 2.9 ) is the imaginary part of a complex- valued 



power series, it is valid in the disc of convergence of this complex- valued series. The 
outermost radius R of this disc is given by 



R = Distance between the point of interest {a^-, 0) 

and the nearest singularity in uj{a); (4.3) 



a sufficient condition for the validity of Eq. (2.9) at the point (ar,ai) is thus 
I ail < R. In most reasonable cases (see Appendix A for a precise formulation of 



these terms), this is also a necessary condition. 

In practice, the requirement that the global topography of the dispersion relation 



w(a) be known will not limit the potential uses of the local formula (2.9): since we 
envisage that our formula should be applied to the in-depth analysis of numerical 
results (for energy-budget analyses and the characterization of mode competition 
and C/A transition curves), this does not result in any loss of relevance for the 
formula (pJl. 



5. Application of the formula to a canonical fluid instability 

The quadratic approximation holds true exactly only for model dispersion relations 



of the kind examined in Eq. (3.7). In more practical systems, the exact series ex- 
pansion ( |2.9[ ) would have to converge quickly for the quadratic approximation to 
be useful. The convergence properties have already been discussed in general terms 
in Sec. El here, we discuss convergence for the particular model (dimensionless) 
velocity field 



Uo{z) = 1-A + 2A{1-Fsinh2^ [zsinh'^l)] }" 



A<0, 



(5.1) 
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where A and N are dimensionless paranieters, and — oo < z < oo. This velocity 
field models the steady wake profile generated by flow past a bluff body, and is 



absolutely unstable (Monkewitz 1988) (this serves as a model test case for our 
purposes; for subsequent work on cylinder wakes, see Pier 2002). The associated 
OS equation possesses three independent parameters, {A,N,Re), where Re is the 
Reynolds number. 

We have solved this equation using the Chebyshev collocation method, equipped 



with basis functions appropriate to the infinite domain (Boyd 2001). We have 
verified the correctness of the method by computing the critical Reynolds number 
for the onset of (convective) instability: Rec{N = 1, A = —1) = 1.8, and Re^N 



2, A = —1) = 1.9, in agreement with the existing literature (Monkewitz[|1988). We 



have further verified our numerical method by solving a similar OS problem, wherein 
artificial confinement is introduced at z = ±7?, such that (j){zLH) = (t)'{±H) = 0. 
This problem is then solved in a standard fashion, using Chebyshev polynomials as 
the basis functions. The value of H is chosen a,s H = H^; here H^ is large in the 
sense that the linear-stability results are the same for H — He and H > H^- The 
only exception is near the imaginary axis, where confinement features occur; these 
features are well separated from the physical saddle points in the dispersion relation 
and do not affect the stability results (Non-physical singularities are identified as 
they move around upon changing the amount of artificial confinement). 



Before applying Eq. (2.9), we investigate the feasibihty of so doing, that is, we 
examine the topography of w(a) and the proximity of any singularities in this func- 
tion to the point of interest (a,-, 0) on the real line. This is done in Fig. [I] The study 
in Fig. [l] is carried out at Re — 100, A = —1.1, and A^ = 2,5. These cases corre- 
spond to absolute instability. In each case, a sequence of large maxima is visible 
along the imaginary axis. We have verified that the sharp maxima along the imag- 
inary axis correspond to poles: successively decreasing the resolution in the scan 
of the complex a-space leads to larger numerical maximum values of (wr,a;i). This 
behaviour is observed in the numerical code involving basis functions appropriate 
to an infinite domain, as well as in the numerical code with artificial confinement. 
By replicating the study using the OS solver with confinement at a large but finite 
value H ~ He, we have verified that these features are due to pseudo- confinement 
in the numerical solver: even with the basis functions appropriate for the infinite 
domain, the numerical solution has artificially interpreted the problem as though 
its geometry were confined. However, the physical saddle point (that is, one not 
due to confinement) is far from the confinement poles on the imaginary axis, and 
hence the latter have no effect on the stability properties of the system. 

The sharp boundaries in the plots for uji- in Fig. IT] indicate mode competition, 
where the most-dangerous mode jumps between distinct eigenvalues of the OS equa- 
tion. However, the ancillary OS study with confinement at a large but finite value 
H = He has demonstrated that the isolated poles along the imaginary axis in 
fact connect with the eigenvalue associated with the temporally most-dangerous 
mode. Thus, the singularity in the most-dangerous mode nearest to the real axis 
lies at ai — —0.12 for both N — 2 and A^ = 5. The physical saddle is located at 
(ar,ai) = (0.79, -1.0) and (0.94,-0.7) for A^ = 2 and Af = 5 respectively, and 



the radius of convergence defined by Eq. (4.3) is therefore R = 1.0 and R = 0.87 



respectively. It would therefore appear that the series (2.9) converges for the N = 2 



case and diverges for the N = 5 case; however, we have found that high-order 

Article submitted to Royal Society 



Temporal and spatio-temporal growth rates: an analytical connection 11 




fi 



-1.5 



0.2 0.4 0-6 O.S 1 1.2 
a 

r 

(a) uJr,N = 2 




0.2 0.4 0.6 S I 12 

(b) a;i,Af = 2 



-0.5r i i 




0.2 0.4 0,6 O.S 1 1.2 
a 

(c) U)r,N = 5 



-0.5 




0-2 0.4 0-6 0.8 I 1.2 
a 

(d) aji,Af = 5 



Figure 1. The global topography of {iu,.,uji) for N = 2 and N — 5. Here Re — 100 and 

A = -1.1. 



truncations of the power series for LUi (and an analogous series for cj^) give good 
approximations to (wj., Wi) outside of this region (Fig. ^. The numerical differen- 
tiation of the Qfr— derivatives for Fig. [2] has been performed by interpolating the 
temporal dispersion relation on to a Chebyshev collocation grid and then perform- 
ing analytic differentiation on the Chebyshev basis functions of this interpolation. 



Possible technical reasons for the apparent success of the power series (2.9) outside 



of its radius of convergence implied by Fig. [2]are discussed in |Appendix A for phys- 
ical applications such as the present one, it suffices to determine on a case-by-case 
basis whether finite truncations of the series for w; are good approximations to the 
underlying function. 

Motivated by the good agreement between the true function LUi(ar,ai) and fi- 
nite truncations of its Taylor series in Fig. [2] we have investigated the temporal 
stability properties of the system (Fig. [3| with a view to applying Eq. (2.9) to 
the spatiotemporal problem. The growth rate has the basic features of a quadratic 
equation in ofr, while the group velocity depends only weakly on the wavenumber 
a^. However, the second and first derivatives of these quantities are non- linear and 
non-constant respectively, and the assumptions of the quadratic approximation do 
not apply. We therefore examine both the quadratic approximation and further. 



higher-order approximations; the latter involve a truncation of Eq. (2.9 1 at cubic 
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Figure 2. Successive truncations of (a) the series for oji (Eq. (2.9l); (b) an analogous series 
for Wr, and a comparison with the true formulae. Here a^ — 0.7729 (i.e. a point close to 
the saddle), Re — 100, A'^ — 5, and A — —1.1. The solid vertical line indicates the radius 
of convergence, while the broken vertical line indicates the ai-location of the saddle point. 





(a) 



(b) 




(c) 




(d) 



Figure 3. Dispersion curves for two test cases with A — —1 and Re — 100. Broken line: 

N =1- sohd line: iV = 8. 



(or higher) order, together with a procedure analogous to Eqs. (3.2)-(3.6) to obtain 
the criterion for the onset of absolute instabihty. 

The true boundary between the convective and absolute regions is calculated 
numerically from a spatio-temporal OS analysis and is identified according to steps 
(i)-(ii) in Sec.fl] The estimated boundary is calculated, according to quadratic and 
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(a) 



(b) 



Figure 4. Boundary between the absolute and convective regions for different Reynolds 
numbers, with comparison to the (a) quadratic, and (b) cubic approximations. Dashed 
curve: Re — 20; solid curve: Re — 100. Circles: approximate curves at Re = 20; crosses: 
approximate curves Re = 100. 



cubic approximations. The quadratic approximation yields a maximum error of 20% 
with respect to the numerically-generated curves. For the cubic approximation, the 
Re = 20 approximate curve and the numerically-generated curve virtually coincide, 
while the maximum error in the Re = 100 curve for the cubic approximation is 
reduced to 7.5% (Fig.[4|. We have verified that going to higher-order truncations in 
Eq. (2.9) for the Re — 100 case leads to near-perfect overlap between these curves 
(a maximum error of 2.5% in a fifth-order truncation, at Re — 100). It should 
further be noted that a second temporally unstable mode exists in the Re = 100 



case: this fails to satisfy the approximate criteria (3.2|-(3.6) for absolute instability 



and, moreover, the numerically-generated dispersion relation for this mode fails to 
exhibit the behaviour necessary for absolute instability. This is in agreement with 
earlier work ^Monkewitz^ (1988) . 

The quadratic approximation can now be used to further analyse the C/A tran- 
sition. The coefficients in a quadratic approximation for the temporal dispersion 
relation uJi{ai) can be related to the temporally most dangerous mode (amjWj^"^), 
such that 



W: 



temp 



2oj. 



temp 



a 

ttm 



a 
am 



quadratic approximation. 



(5.2) 



Moreover, we have found in our numerical calculations that the ar-location of the 
unstable saddle point in the complex dispersion relation almost coincides with the 
most dangerous mode, such that (using the notation in Sec. ^b)) a* ~ am. The 
criterion (3.6) in Sec.[3jb) describing the onset of absolute instability then reduces 
to 



2w 



tcinp 



r(am). 



(5.3) 



Both sides of the condition (5.3) are well represented by expressions of the form 
A{N, Re)+B{N, Re) A over the part of parameter space considered throughout this 
section. In Fig. [5] the dependencies of these coefficients on TV and Re are shown. 
The most significant of these is the TV— dependency of the coefficient B for the 
group velocity: increasing N makes the group velocity approach unity rapidly from 
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Figure 5. Dependency on A'^ of coefficients in fits linear in A of the left-hand side of Eq. ( 3.6 1 
(a) and the group velocity (b) at Re = 100 (solid lines) and Re = 20 (dashed lines); plusses 
and crosses represent A and B, respectively. In (a), the sign of both coefficients has been 
reversed. 



below. If we ignore for the time being the other dependencies on N and use for 
the left-hand side of Eq. (5.3 1 simple A^— averaged values (A = oq ~ —0.05 and 
B = ai~ —0.8 for Re = 100), and for the group velocity the curve fit 



C(,(a,„) « 1 -f- feoexp(-6iA^)A 



(5.4) 



(with 6o ~ 0.6 and hi ~ 0.5 for Re = 100), we obtain an approximate expression 
for the C/A transition based on the quadratic approximation: 



A 



CA 



flo - 1 



boexp{—biN) — oi 



(5.5) 



We have found this to represent the full results obtained with the quadratic ap- 
proximation within the uncertainty of the value of the coefficients. Overall then, the 
main trend is that in making A more negative, the group velocity (the right-hand 
side of Eq. (5.3)) decreases, whereas the left-hand side of Eq. (5.3) increases; these 
two effects promote absolute instead of convective instability. Thus, the steep drop 
in the curves at large N in Fig. |4] is associated with the fact that increasing N 
increases the group velocity. 

Regarding the dependencies on Re, we first note that at low values of N, the left- 
hand side of Eq. (5.3) is small at Re = 20 compared with Re = 100, whereas right- 



hand side is large; these changes both favour convective instability. At intermediate- 
to-large values of N, both terms decrease when the value of Re is decreased, with 
the left-hand side more so than the right-hand side. Hence, the C/A transition 
in Fig. |4] occurs at a more negative value of A for the lowest value of Re due to 
effects on both sides in Eq. (5.3). These effects can be accounted for in the above 



through the introduction of explicit i?e— dependencies in the coefficients aq, ai, 6oj ^i 
(even if only representing the i?e— dependencies of their A^— averaged values). These 
dependencies are already visible in Fig. [5l although their parametrization is beyond 
the scope of this paper. 

These somewhat crude approximations can be refined by representing the de- 
pendency of the left-hand side coefficients A and B on N, the causes of which 
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Figure 6. The dependencies on A^ and A of (a) — d^cjj'^^^/dar ; (b) o;''^™^, at the most dan- 
gerous temporal mode for Re = 100. Results are shown for TV = 1, 1.25, 1.5, 2, 3, 4, 5, 6, 7, 8 
from top to bottom in (a), and from bottom to top in (b). 



we first briefly investigate here. At values of N near unity, the left-hand side of 
Eq. (5.3 1 first increases strongly with N before decreasing beyond A^ « 3, such that 



the region of absolute instability in parameter space is reduced further at large 
values of N. In fact, the results in Fig. l5Fa) tell two different stories, for values 
below and above iV = 3. Referring back to Eq. (3.6 1, we recall that the left-hand 



side is the product of the maximum temporal growth rate and minus the curvature 
of the dispersion relation. Results for these two components are shown in Fig. |6] 
We conclude from Fig. 6 that u!^°"^^ {a^^) increases uniformly with increasing N 



consistent with the increase in the entire left-hand side of Eq. (3.61 at low values 
of N, and that the resulting promotion of absolute instability is only reversed by a 
decrease in —d'^uj^'^'^^/da^ at large values of N; at low values of N, this curvature 



term hardly depends on N. It is possible to account for these trends in Eq. (5.51 



by using bilinear fits for 0;;,^^ and a,^ (as well as Eq. (5.4) for the group velocity) 



temp 



«a + alN + b^A + blNA, 
a" + alN + 6° A + 5,1, A^A. 



(5.6) 



Substitution of Eq. (5.61 into Eq. ( 5.3 1 then yields a quadratic equation for A at the 



C/A transition. This whole approach can be extended by lifting the approximation 
a* ~ am, or by going to higher-order truncations of the dispersion relation; however, 
these complications obscure the understanding achieved here using simple linear and 
bilinear fitting, and such extensions are not pursued further here. 



6. Discussion 

We have presented an analytical connection between spatio-temporal and temporal 
growth rates in a linear stability analysis. Specifically, we have derived approximate 
criteria for a transition between convective and absolute instability, explicitly in 
terms of the temporal problem, to various levels of approximation. The simplest 
approximation is a quadratic one where the location of the saddle a* is assumed 
to coincide with the temporally most-dangerous mode; then, the criterion for the 
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onset of absolute instability is given by Eq. ( 5.3 ). At the next level of approximation 
(still quadratic), one determines a* via Eq. (3.51 Although this approach can be 
used directly to determine C/A transitions, we anticipate that the main use of 
the theory will be in analyzing results obtained with an efficient computational 
algorithm for the full OS problem. We imagine that one would perform such a 
numerical calculation, diagnose absolute instability (using the saddle-point/pinch- 
point criteria) , and then characterize the instability in detail. Our formula can help 
in this characterization in the following ways: 

1. In parametric studies, the dependency of C/A transition curves on model pa- 
rameters can be difficult to analyse, especially if there are many governing 
parameters. In Sec. [5] we have seen an example wherein simple approximations 
based on Eq. ( |2.9[ ) , together with insight into the trends exhibited by the dis- 
persion relation in the purely temporal problem, yield a detailed explanation 
of the behaviour of the C/A transition curves. 

2. Temporal mode competition is observed in various systems and, as pointed 
out in the Introduction, disentangling different spatio-temporal modes poses a 
formidable problem. The present approach then provides a convenient aid: for 
each temporal branch, an approximation of the corresponding spatio-temporal 



growth rate in the complex a plane results directly from Eq. ( 2.9 ) . An example 
is given in Fig. [2] Although only a cross-section for a single value of a-^ is 
shown in Fig. [2] this is readily extended to entire sections of the complex a 
plane. Indeed, in computing the C/A transition in truncations higher than 
quadratic, we have found this approach more straightforward rather than the 
root- finding approach taken for the relatively simple quadratic case (Sec.p|. 

3. It follows from the previous point that features in contour plots for a spatio- 
temporal growth rate (such as saddle points) can be related directly to a 
specific temporal mode with the present approach. The advantage of this is 
that such features can then be associated with the physical mechanism that 
governs that temporal mode. 

The appl ication of the approach p roposed here to real flows has recently been com- 



pleted in O Naraigh et al. (2012), where we have studied absolute and convective 
instability in two-phase gas/liquid flows. There, however, the focus was on the appli- 
cation of the theory in a detailed parameter study and a derivation of the quadratic 
approximation was not included; in the present work, the focus is on the mathemat- 



ical derivation of Eq. (2.9), discussion of the radius of convergence of the formula. 



and the investigation of the accuracy of finite truncations of the equation (2.9). 
Such analysis is important before the theory can be used with confidence in further 
applications. 

Throughout this paper we have considered two-dimensional disturbances, taking 
w as a complex- valued function of a single complex variable (the complex wavenum- 
ber). However, the approach taken could be extended to three-dimensional spatio- 
temporal disturbances, via a double- Taylor-series expansion of the dispersion rela- 
tion in the streamwise and spanwise wavenumbers. For cases wherein the spanwise 
wavenumber is real, the analysis presented herein requires no modification, as the 
spanwise wavenumber would then only enter as an additional parameter. 
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We conclude by discussing briefly a further application of this method to weakly 
globally unstable modes. In nonparallel flows, the properties of the system vary in 
the streamwise direction, and a standard normal-mode decomposition is no longer 
possible. The generic linear-stability problem to be solved then reads 

dtp r f ^ ^ 
dt V 9x ' dz 



v-, 



where £ is a linear operator with coefficients that depend on the streamwise direc- 
tion; the normal-mode decomposition d/dx — ^ ia is no longer possible. We make 
trial substitution -0 — >■ i^{x, z)e~'""^* and obtain the global eigenvalue problem: 

d_ d_ 
dx' dz ^ 



'iiUQip{x, z) = C 



ip{x,z) 



The system is globally unstable if there exists a global mode for which S(wg) > 0. 
We introduce a parameter e = X/L, where A is the typical length scale of the 
disturbance ijj, and L is a measure of the system's extent in the streamwise direction. 
If £ <C 1, and if the coefhcients of the linear operator C depend only weakly on 
space, through the combination X := ex, then a WKB approximation (Huerre and 



Monkewitz 1990 Chomaz et al. 1991 Schmid and Henningson 2001 ) can be used 



to demonstrate that the leading-order approximation to the most-dangerous global 
mode is 

WG=wo(^s) + C'(e), 

where ujq{Xs) is the saddle-point frequency calculated in the usual manner from 
a purely local stability analysis, at the fixed parameter value Xg, and Xg is the 
location of the saddle point of the complex- valued function u!q{X) in the complex 
X-plane: 

dujQ 



dX 



(Xs) - 0. 



The flow is therefore globally unstable if luq{Xs) > 0. Thus, we are interested in 
a prescribed complex- valued function of a complex variable, ujq{X), and its saddle 
point(s). If uJo{X) is holomorphic, and if ujo{X) is known along the real line, then 
the arguments of Sec. [2] apply also to the function ujo{X), and the following formula 
pertains: 



Woi(^r, ^i) = WOi(^r, Xi=0)+Yl 



(-1)" d^n-hl 



(2n + l)!dX2"+ 



j[uJoAXr,Xi = 0)]X( 



2n-|-l 



E 



(-1) 



n+l 



j2n+2 



(2n + 2)!dx2"+2 



[uJoi{X,,Xi = 0)]X( 



2n+2 



(6.1) 



Furthermore, if uJo{X) is quadratic in X, then the truncation in Sec. 
and the condition for the onset of global instability reads 



c D 



r^2 



d^LOQi iX,, Xj = 0) 
dX} 



woi(^;,^i-o) 



xt 



da;or(^r,^i ==0) 



dX, 



applies. 



(6.2) 



x: 



where X* satisfies a root-finding condition analogous to Eq. (3.5). It is anticipated 



that this application of our formula to global stability analysis, together with the 
techniques (l)-(3) for analyzing the results of local stability calculations, will find 
application in stability problems beyond those mentioned in this work. 
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Appendix A. Radius of convergence and optimal truncation 

In this Appendix we demonstrate that the usual criterion for the convergence of 
a complex power series applies to the central result of this paper, namely the se- 



ries (2.9). It is tempting to conclude that the convergence properties of the re- 
sult (2.9 1 do not depend on the global topography of the function w(a), and that 
a local analysis (similar to that used in Taylor's theorem in real analysis) applies. 
We demonstrate here that this conclusion is incorrect. We also discuss the notion 
of an optimal truncation for a divergent power series, such that a finite truncation 
of a divergent power series can be used to approximate the generating function 
of the power series. We demonstrate the applicability of this method to the series 
examined in Sec. [5] thus explaining the apparent 'over prediction' obtained in that 
section. 

The starting-point for the convergence analysis is Taylor's theorem applied to 



the purely real Taylor expansion ( 2.9 ). We introduce some notation (see also Fig.^: 



• 



• 



The point O = (aro, 0): The centre of the power series based on the real axis. 

(aiOiCCio): The point of interest at which the power series is to be evaluated. 
For our purposes, aio will be negative. 

The point C = (ctrOj otic)'- A point used to compute the remainder in Taylor's 
theorem; here a^c € [0, aio]- 



• The point X = a^: The singularity closest to C. 

We shall assume that the derivatives ((9''wi/9af) exist for all orders k S 
{0, 1, 2, • • • }, on the interval (0, aio) (the brackets (•, •) here denote an open inter- 
val, and not a point in the complex plane). Taylor's theorem guarantees that in 



approximating aji(Q;r,cki) by a finite truncation based on Eq. (2.9), we incur an 
error 

1 9"+iwi 



Rn 



(n+1)! aar+i 
Using Cauchy's theorem, we have 



n+l 
iO ' 



aro+\ aice[0,aio]. (Al) 



Qn+l^ 



da 



n+l 



d"+lc 






(n + l)! [ uj{a) 



(Q,o,aic) 2^^^ Jl [" - ("rO + ittic)]" 



da, (A 2) 



where 7 is a circle centred at (aro, aic); of radius R. Here we have switched from 
a partial derivative to a standard complex-valued derivative, because the function 
uj{a-[, a{) is analytic away from singularities, and the limits calculated in respect of 
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Figure 7. Figure showing the centre of the power series, the singularity, and the point of 
interest at which the series is to be evaluated. 

the derivative are independent of approach. The radius R is yet to be determined. 
Also, 



a"+iwi 



dal 



n+l 



(arO:"ic) 



< 



gn+1^ 



da 



n+l 



(oirO-aic) 



(n+l)! / 

< max|aj(Q;) 

R V 7 



(A3) 



Combining Eqs. (All and (A3|, we have 



\Rn\ < ( max \uj\ 
V 7 



"io 



R 



n+l 



To make |i?„| — >■ as n — )■ cxd, we would hke to take R > |aio|- However, we cannot 
take R as arbitrary, we must have R < \CX\, where \CX\ is the distance between 
C and the closest singularity thereto. Because we have no a priori knowledge of 
the location of C, we must take account of the worst-case scenario, where C = O. 
Then, for the series to converge, we must have 



R< \0X\, 



(A 4) 



that is, in order for the power series centred at O to converge, the point of interest 
(arOiCtio) must be inside a disc of radius R, where R is the distance between the 
power-series centre and the nearest singularity. In other words, the standard crite- 
rion for complex power series applies: the power series is valid provided the point 
of interest lies within the disc of convergence. 

Equation (A4) is a sufficient condition for the Taylor series for uji{a^,ai) to 



be valid. If we extend this Taylor series using the Cauchy-Riemann conditions 
and construct a Taylor series in ai for Wr(Q;r, ai), the standard results of complex- 



variable theory apply, and Eq. ( A 4 1 is both sufficient and necessary in order for 



the Taylor-series pair for (a;i.,a;i) to converge. On the other hand, it is possible 
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for the Taylor series of Wi to converge beyond the radius of convergence, while at 
the same time the Taylor series for Wr diverges. However, this rather pathological 
scenario is unlikely to occur in the applications contained in this paper: for example. 
Fig. [2] shows that the finite truncations to the power series for both uj{ and uj^ are 
good approximations to the underlying dispersion relation outside of the radius of 
convergence of the power series. Our final discussion therefore concerns the notion 
of optimal truncation for divergent power series. 

Consider a real-valued function f{x) that admits a power-series expansion, cen- 
tred at the origin: 

/p(:r) = Vc„x", c„ = ^^-^. (A5) 

Assume moreover that the power series fp{x) converges in a finite interval |a;| < R, 
and that fp{x) — f{x) on the same interval. If fp{x) remains finite at a; = i?. 



then it is possible to use a finite truncation of Eq. (A5) to approximate f{x) in 
the region \x\ > R. In other words, for a given x such that |a;| > R, there exists a 
positive integer K such that the difference 

\ n=0 / 

is minimized. As an example of this sort, consider the function f{x) — (1+a;) log(l-|- 
x) — X. On the interval \x\ < 1, this function can be represented as a power series: 

°° n+l 

/(x)^^(-l)«+i-^--. (A 6) 

^-^ nin + 1 

n—l ^ ' 

It is straightforward to show that the power series evaluated at x = 1 converges 
(the limit comparison test applies). Thus, we extend finite truncations of the power 
series beyond x = 1 to approximate the underlying function ]{x). In doing so, 
we wish to minimize the difference bKix). Using elementary facts about geometric 
progressions, this difference can be re- written as 



^k{x) := I / Ax \ Ax" — Y^^' ^^^^ 




2 



The minimum iC-value satisfies 



''&w=o, rd.'rd."'-"'"'°y"'.o. (A 7) 



dK 



JQ 



For a; > 1, Eq. ( |A7[ ) has a solution with positive K. We have plotted the optimal 
value of K in Fig. I8[a) . The difference between the optimal truncation of the diver- 
gent power series and /(a;) is shown in Fig. IsFb). There is little difference between 
the optimal truncation and a uniform truncation with K = 20. The relative error 
between the uniform truncation and the true function /(x) is shown in Fig. ^c): 
this never exceeds 2% on the interval a; G [1, 1.1]. 

Finally, in Fig. |9]we examine the coefficients of the Taylor series for the model 
dispersion relation developed previously (Sec. I5|. We consider the case Re = 100, 
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(a) 



(b) 




(c) 



Figure 8. Analysis of finite truncations of the divergent power series (A6l. (a) 
The optimal truncation value K; (b) The difference between the true function 
f{x) — (1 + a;) log(l + x) — X and the optimal truncation; (c) the relative error in making 
a uniform truncation. 



N = 5, and A = —1.1. Panel (a) refers to ai- ~ 2.1978, and the singularities on 
the imaginary axis he far from the point of interest [a* , 0) (here a* refers to the 
real coordinate of the saddle point). The first eight coefficients in the power series 
exhibit exponential decay. On the other hand, in panel (b) (for which a^ — 0.7739), 
the coefficients decay algebraically. This case corresponds to a ar ~ a* (i.e. close to 
the saddle), for which the radius of convergence is i? = 0.87; in the main part of the 
paper we extended finite truncations of the series beyond this point, i.e. |ai| > R. 
That the first eight coefRcients of the power series decay algebraically suggests that 
the continuation described herein is appropriate. Thus, it is possible (albeit with 
great caution) to use a finite truncation of a divergent power series to approximate 
the generating function of the power series. 
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Figure 9. Coefficients of the Taylor series for (2.91 {A„) and for the analogous series in tJr 
(Bn), for the model discussed in Seep] Here, Re = 100, and N = 5. Case (a) corresponds 
to «! = 2.1978; case (b) corresponds to a,- = 0.7739 (i.e. a point close to the saddle, for 
which the radius of convergence is small {R = 0.87)). 
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